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SUMMARY 

The trajectories of droplets in the air flowing past an NACA 65A004 
airfoil at an angle of attack of 8° were determined. The amount of 
water in droplet form impinging on the airfoil, the area of droplet im- 
pingement, and the rate of droplet impingement per unit area on the air- 
foil surface were calculated from the trajectories and presented to cover 
a large range of flight and atmospheric conditions. These impingement 
characteristics are compared Briefly with those previously reported for 
the same airfoil at an angle of attack of 4°. 


INTRODUCTION 

The data presented herein are a continuation of the study reported 
in reference 1 on the impingement of cloud droplets on a low-drag, thin 
airfoil. The airfoil studied in both the reference cited and in this 
report is a 4-percent-thick symmetrical NACA 65A004 airfoil. In refer- 
ence 1 the impingement characteristics of the airfoil were reported with 
the airfoil set at an angle of attack of 4 ; whereas, the data herein 
apply for an angle of attack of 8°. The data calculated for 8° angle of 
attack, along with those data for 4° presented in reference 1, permit the 
evaluation of the impingement characteristics for circling, landing, and 
some types of flight plans for pursuit or fighter aircraft. 

The trajectories of atmospheric water droplets about an NACA 65A004 
airfoil at 8° angle of attack at subsonic velocities were calculated 
with the aid of a differential analyzer at the NACA Lewis laboratory. 

From the computed trajectories, the rate, distribution, and surface ex- 
tent of impinging water were obtained and summarized in this report. 


SYMBOLS 

The following symbols are used in this report: 
d droplet diameter, microns (micron = 3.28X10"® ft) 
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L 


Re, 


U 

u 

W 


W_ 


m 


W Q 


w 


x>y 


P 


inertia parameter, 1.704X10 -12 dimensionless (the density of 

water, 1.94 slugs/cu ft, is included in the constant) 

airfoil chord length y ft 

—6 

free-stream, Reynolds number with respect to droplet, 4.813X10 
dp a U 

, dimensionless 

P 

distance on surface of airfoil measured from leading-edge chord 
point, ratio to chord length 

flight speed, mph 

local air velocity, ratio to free-stream velocity- 

rate of water impingement per unit span of airfoil, lb/(hr)(ft span) 

rate of total water impingement per unit span of airfoil, lb/(hr) 

(ft span) 

local rate of water impingement, lb/(hr){sq ft) 
liquid-water content in cloud, g/cu m 
rectangular coordinates , ratio to chord length 


local impingement efficiency, 


% 
dS ’ 


dimens ionle s s 


(jl viscosity of air, slugs/(ft) (sec) 
p density of air, slugs/cu ft 

cl 

Subscripts : 

l lower airfoil surface 

s airfoil surface 

u upper airfoil surface 


0 


free stream 



NACA TN 3155 


3 


RESULTS AND DISCUSSION 

In order to obtain the extent of impingement and the rate of droplet 
impingement per unit area on the airfoil, the cloud droplet trajectories 
with respect to the airfoil were determined. The method for calculating 
the droplet trajectories is described in reference 2. A solution of the 
differential equations that describe the droplet motion was obtained with 
the use of the mechanical analog (described in ref. 3) based on the prin- 
ciple of a differential analyzer. The air-flow field around the airfoil 
was obtained by the vortex substitution method described in reference 2, 
except that, for the NACA 65A004 airfoil, the velocities at the surface 
of the airfoil were calculated by the method described in reference 4, 
whereas the surface velocities on the 65-^212 and 65-^208 airfoils dis- 
cussed in reference 2 were obtained from wind-tunnel measurements of the 
pressure coefficients. The values of the surface velocities for the 
65A004 airfoil were calculated by the Douglas Aircraft Corporation for 
the Lewis laboratory (see fig. l) . Although the droplet trajectories 
were calculated for an incompressible flow field, the results of the 
calculations can be applied up to the flight critical Mach number (ref. 5). 

The geometric chord line of the airfoil is oriented at an angle of 
8° with the x-axis of the rectangular coordinate system, and the leading 
edge is placed at the origin of the coordinates, as shown in figure 2. 

The airfoil orientation presented in references 1 and 2 is retained here- 
in, except for the magnitude of the angle of attack. At an infinite 
distance ahead of the airfoil, the uniform air flow carrying the cloud 
droplets is assumed to be approaching the airfoil from the negative x- 
direction and parallel to the x-axis. All distances are dimensionless, 
because they are ratios of the respective actual distance to the airfoil 
chord length L . 


Rate of Water Interception 

The rate of total water interception, in pounds per hour per foot of 
wing span, is determined by the tangent droplet trajectories (fig. 2), 
by the speed of the aircraft, and by the liquid-water content in the 
cloud. The flight speed and size of the airfoil, as well as the droplet 
size in the cloud, are the principal variables that affect the spacing 
between the two tangent trajectories. The amount of water that strikes 
the airfoil is proportional to the spacing y Q ^ - y Q j, and the rate of 

total water interception per unit span of the airfoil on that portion of 
the airfoil surface bounded by the upper and lower tangent trajectories 
can be calculated from the relation 

W m = 0.33UwL(y 0jU - y Q>l ) 


( 1 ) 
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The values of y - y are given in figure 3 in terms of the rec ip- 

0, u 0,2 

rocal of the inertia parameter l/K and the free-stream Reynolds number 
Re 0 * The inertia parameter K is a measure of the droplet size, the 

flight speed and size of the airfoil, and the viscosity of the air 
through the relation 


2 

K = 1.704X10' 12 (2) 

The density of water and the acceleration of gravity, which are expressed 
as part of the conversion factor, are 62.4 pounds per cubic foot and 
32.17 feet per second per second, respectively. The free-stream Reynolds 
number is defined with respect to the droplet as 

Re n = 4.813X10' 6 (3) 

A graphical procedure for determining values of the dimensionless param- 
eters K and ReQ in terms of airplane speed, chord length, altitude, 

and droplet size is presented in appendix B of reference 2. 

The variation of rate of water interception with airfoil speed is 
summarized for an altitude of 20,000 feet in figure 4, in which the ordi- 
nate W m /w is the total rate of water impingement per foot span of air- 
foil per unit liquid -water content (g/cu m) in the cloud. Several chord 
lengths ranging in value from 2 feet to 20 feet are considered. The val- 
ues in figure 4 are for flight through clouds composed of uniform drop- 
lets 15, 20, 30, and 40 microns in diameter. The values of W /w. given 

in figure 4 are based on the most probable icing temperature as a func- 
tion of altitude presented in figure 15 of reference 2. (The most prob- 
able icing temperature was obtained from approximately 300 icing obser- 
vations in flights.) As shown in reference 2, a change in altitude of 
10,000 feet will change the rate of water impingement by approximately 7 
percent. The droplet size and the liquid-water content of clouds are 
seldom known with sufficient accuracy (ref. 3) to permit the rate of 
water impingement to be calculated within 10 percent; therefore, within 
practical limits of application, the results of figure 4 can be used over 
a wide range of altitudes (approx. ±10,000 ft, see ref. 2). 

The effect of wing taper can also be obtained from figure 4, provid- 
ed that for each section of span considered the taper is small enough 
that two-dimensional flow over the section is approximated, as is men- 
tioned in reference 2. 
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Extent of Impingement 

The limit of impingement is determined by the point of tangency on 
the airfoil surface of the two tangent trajectories. The rearward limits 
of impingement on the upper surface are shown in figure 5(a), and on the 
lower surface in figure 5(b) . The distances S u and Sj are measured 

on the airfoil surface from the point of intersection of the geometric 
chord line with the leading edge (fig. 2) in terms of the chord length. 

The limits of impingement are given in figure 5 in terms of the recipro- 
cal of the inertia parameter and free-stream Reynolds number. 

At an angle of attack of 8° the extent of impingement on the upper 
surface is always less than 2 percent of chord (S u = 0.02), which is the 

largest value possible for the extreme conditions where l/K = 0. The 
value of S u for l/K = 0 was obtained from simple geometric relations. 

The remaining values of presented in figure 5(a) were obtained from 

calculated trajectories . These calculated values may be quite inaccu- 
rate, because of the sharp-edged shape of the airfoil at the leading edge. 
The accuracy of determination of S u is unimportant in an application 

to airfoil icing calculations, because all values of S u must be less 

than 0.02. The estimated accuracy for the values given in figure 5(a) 
is +20 percent . 

Nearly all the water impinging on the airfoil impinges on the lower 
surface between the leading-edge chord point and the lower limit given 
in figure 5(b) . As was discussed in reference 1 for the NACA 65A004 air- 
foil at 4° angle of attack, the tangent trajectories approach the lower 
surface of this airfoil very gradually. This very gradual approach of 
the lower-surface tangent trajectories leads to uncertainties as to the 
location of the tangent point, and the values of Sj are therefore sub- 
ject to individual interpretation. The uncertainties in the location of 
the tangent point can range from ±4-percent chord at l/K hi to +7 
percent at l/K — 3 and back to +4 percent at l/K - 100. For the ex- 
treme case of l/K = 0, the impingement extends to the trailing edge of 
the airfoil section. The lower-surface limits are summarized in figure 
6 for the same speeds, chord lengths, droplet sizes, and altitude given 
in figure 4. 


Impingement Distribution on Surface 

Trajectory starting ordinate as function of point of impact . - The 
manner in which water is distributed on the surface of an airfoil can be 
obtained if the starting point of a droplet trajectory is known with 
respect to the point of impingement on the surface. The starting 
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ordinate 


'0 


at infinity of any impinging trajectory, including trajec- 


tories bounded by the upper and lower tangent trajectories (fig. l) , can 
be found in figure 7 with respect to the point of impingement on the sur- 
face . The values for the starting and ending positions of the trajec- 
tories are shown in figure 7 for four values of free-stream Reynolds num- 
ber. For each value of ReQ, curves for several values of l/K are 


given . 


The amount of water impinging between any two given points on the 
airfoil surface may be found by applying the results of figure 7 in the 
relation 

W = 0.33UwL(y 0 x - y Q ^) (4) 

The values of y^ - Jq ^ obtained from the end points of each 

curve in figure 7 are the same as the values given in figure 3. The 
value of y Q ^ for l/K = 0 (not shown in fig. 7) is -0.139 at 

S 7 = 1.00229 (airfoil trailing edge) for all values of Re Q . 

L 


Local rate of droplet impingement . - The local rate of droplet im- 
pingement per unit area of airfoil surface can be determined from the 
expression 


W 


P 


<a-y r 


0 . 33Uw- 


dS 


0 . 33UwP 


(5) 


which is related to equation (4), with proper consideration for the fact 
that y 0 and S are based on the wing chord L. The values of the 

local impingement efficiency P as a function of the airfoil distance 
S are given in figure 8. These values were obtained from the slopes of 
the curves in figure 7 . 

As is discussed in reference 1, the values of p (fig- 8) are very 
sensitive to the shape of the y Q against S curves (fig. 7). Because 

of the geometry of the sharp-nosed NACA 65A004 airfoil and the manner in 
which the trajectories approach the airfoil surface, small errors in the 
calculated trajectories result in considerable error in the slopes of 
the curves of figure 7. The possible error in the values of p, due to 
the computational procedure, for surface positions other than near the 
stagnation point, is estimated in reference 1 to be somewhat less than 
±10 percent for the values reported therein. Because of improved tech- 
niques in the computational procedure, the values of P given in figure 
8 herein are in error by somewhat less than ±2 percent for surface posi- 
tions other than within 1 percent of the surface distance where the peak 
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values of 8 occur. Since the total water impinging is directly relat- 
ed to 



( 6 ) 


a check on the computational accuracy of the values of 8 in figure 8 
was also obtained by comparing the area under each 8 curve with the 
values of y n - y n given in figure 3. The area values checked with- 

U,U U y if 

in +1.5 percent of the corresponding values for total water interception. 


The accuracy in the maximum value of 8 is also improved in this 
report as compared with that reported in reference 1. The possible error 
was estimated to be +25 percent for the maximum values of 8 reported 
in the reference cited; whereas, the possible error in the maximum val- 
ues of 8 given in figure 8 is estimated to be less than +12 percent. 

As was discussed in reference 2, this possible error is not considered 
very serious, because only a small portion of the total water impinging 
on the airfoil is involved in the error. 


Comparison of Impingement at 8° with Impingement at 4° Angle of Attack 

Rate of total water interception . - For all values of the reciprocal 
of inertia parameter and free-stream Reynolds number, the rate of total 
water interception is greater at an angle of attack of 8° than at an 
angle of attack of 4°. This comparison can be made between figure 3 of 
this report and figure 3 of reference 1. The comparison is summarized in 
the following table for conditions established at 300 miles per hour, a 
chord length of 9.4 feet, and an altitude of 10,000 feet: 


Reciprocal 
of inertia 
parameter, 
l/K 

Droplet 

diameter, 

a, 

microns 

Free-stream 
Reynolds 
number , 

Re o 

yo,u - 

yo,z 

Ratio of rate of 
total water at 
8° to rate at 4° 

Angle of 

attack 

8° 

4 o 

1 

80 

594 

0.099 

0.053 

1.87 

10 

25 

190 

.043 

.021 

2.05 

100 

8 

59 

.017 

.005 

3.40 


The value of free-stream Reynolds number tabulated is the value obtained 
from equation (3) for the physical conditions established for this 
comparison. 
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Extent of impingement. - Although the limit of impingement on the 
upper surface is farther rearward for an angle of attack of 4^ than for 
8°, the total extent is very small in "both cases. At 4° angle of attack 
S u is always less than 0.02 for values of l/K>l. At 8° angle of 

attack S u is always less than 0.01 for values of l/K>l. 

The impingement on the lower surface extends hack to the trailing 
edge at both angles of attack for l/K = 0. The lower-surface limit is 
summarized in the following table for the same flight and atmospheric 
conditions given in the preceding section on rate of total water inter- 
ception: 


Reciprocal 
of inertia 
parameter, 
l/K 

Droplet 
diameter , 

Free- stream 
Reynolds 
number , 

Limit of impingement 
on lower surface, Sj 

microns 

Angle of attack 

u 

8° 

4° 

1 

80 

594 

0.74 

0.53 

10 

25 

•190 

.34 

.17 

100 

8 

59 

.21 

.03 


Local rate of droplet impingement . - At a given point on the lower 
surface, the local rate of droplet impingement p is greater at an angle 
of attack of 8° than at 4°. For both 8° and 4°, the maximum rate of 
local impingement occurs between 0 and S = 0.01 on the lower surface. 
The peak values of P are higher at 8° than at 4°. At 8° angle of 
attack the peak values are between 0.8 and 1.0; whereas, at 4 the peak 
values are reported as low as 0.3 and none higher than 0.56. These large 
differences cannot all be accounted for by the possible estimated error 
in determining the peak values . 


CONCLUDING REMARKS 

The data presented herein apply directly to flights in clouds com- 
posed of droplets that are all uniform in size and to nonswept wings of 
high aspect ratio. A detailed procedure for weighting the impingement 
of droplets for flights in nonuniform clouds is presented in reference 3. 
A method for extending the impingement calculations for nonswept wings 
to swept wings is presented in reference 6. As is discussed in refer- 
ence 5, the impingement results should be applicable for most engineering 
uses throughout the subsonic region, because the subsonic compressibility 
of air does not affect the droplet trajectories appreciably. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, May 20, 1954 
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Chordwlse station, ratio to chord length 

Figure 1. - Velocities on surface of 65A004 airfoil. Angle of attack, 8°; incompressible flow field. 





Figure 3. - Difference between upper and lower trajectories at free-stream conditions for 65A004 
foil. Angle of attack, 8°. 

















Rate of water impingement per unit liquid-water content, W^w, (lb/(hr)(ft span))/(g/cu 
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Limit of impingement on lower airfoil surface, ratio to chord length 
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Figure 6. - Limit of Impingement along lower surface of 65A004 airfoil. Altitude, 20,000 feet 
angle of attack, 8°; most probable icing temperature, -11° F. 





NACA TN 3155 


19 



Figure 6. - Concluded. Limit of Impingement along lower surface of 65A004 airfoil. Altitude, 20,000 feet; 
angle of attack, 8°; most probable icing temperature, -11° F. 
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Figure 7. - Continued. Trajectory starting ordinates as function of point of impingement on surface of 65A004 airfoil. Angle 
of attack, Q°. 
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Distance on airfoil surface, S, ratio to chord length 



■Reciprocal of 
inertia parameter 



Figure 7. - Concluded. Trajectory starting ordinates as function of point of impingement on surface of 65A004 airfoil 
of attack, 8°. 












(c) Free-stream Reynolds number, 256. 

Figure 8. - Continued. Local Impingement efficiency of NACA 65A004 airfoil. Angle of attack, 8 C 




Figure 8. - Concluded. Local impingement efficiency of NACA 65A004 airfoil. Angle of attack, 8 



